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A systematic investigation of the microscopic conditions stabilizing itinerant ferromagnetism 
of correlated electrons in a single-band model is presented. Quantitative results are obtained by 
quantum Monte Carlo simulations for a model with Hubbard interaction U and direct Heisenberg 
exchange interaction F within the dynamical mean-field theory. Special emphasis is placed on the 
investigation of (i) the distribution of spectral weight in the density of states, (ii) the importance 
of genuine correlations, and (iii) the significance of the direct exchange, for the stability of itinerant 
ferromagnetism at finite temperatures. We find that already a moderately strong peak in the density 
of states near the band edge suffices to stabilize ferromagnetism at intermediate U -values in a broad 
range of electron densities n. Correlation effects prove to be essential: Slater-Hartree-Fock results 
for the transition temperature are both qualitatively and quantitatively incorrect. The nearest- 
neighbor Heisenberg exchange does not, in general, play a decisive role. Detailed results for the 
magnetic phase diagram as a function of U, F, n, temperature T, and the asymmetry of the density 
of states are presented and discussed. 
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I. INTRODUCTION 

In contrast to conventional superconductivity and an- 
tiferromagnetism, metallic ferromagnetism is in general 
an intermediate or strong coupling phenomenon. Since 
there do not exist systematic investigation schemes to 
solve such types of problems the stability of metallic fer- 
romagnetism is still not sufficiently understood. This 
is true even within the simplest electronic correlation 
model, the one-band Hubbard model,El in spite of signifi- 
cant progress made recently. The Hubbard interaction is 
very unspecific, i.e., does not depend on the lattice struc- 
ture or dimension. Hence the dispersion, and thereby 
the shape of the density of states (DOS), is of consider- 
able importance for the stability of ferromagnetism.-. This 
was recognized already by Gutzwiller,El Hubbard, and 
KanamorH in their initial work on the Hubbard model. 
However, the approximations used in the early days of 
many-body theory were usually not reliable enough to 
provide definite conclusions. An exception are the exact 
results by Nagaokafl on the stability of ferromagnetism 
at U = oo in the case of one electron above or below half 
filling. They show an important lattice sensitivity but, 
unfortunately, are not applicable in the thermodynamic 
limit. 

Over the years the stability of metallic ferromagnetism 
has turned out to be a particularly difficult many-body 
problem whose explanation requires subtle nonperturba- 
tive techniques. There has baep^ an . u psurge of inter- 
est in this topic most recentlyBil3't§ll3. These investi- 
gations confirm that ferromagnetism is favored in sys- 
tems with (i) frustrated lattices (which suppress antifer- 
romagnetism) and (ii) high spectral weight near the band 
edge closest to the Fermi energy (which improve the ki- 
netic energy of the polarized electrons). Taken together, 
these properties imply a strongly asymmetric DOS of the 
electrons. Ferromagnetism on bipartite lattices having a 



symmetric DOS may still be possible, but seems to re- 
quire very large values of C/.E3 With the exception of 
Refs. [l3|,[l5| all previous calculations refer to the ground 
state. It is therefore of interest to obtain an answer to the 
question: How does the distribution of spectral weight in 
the DOS influence the stability regime of ferromagnetism 
at finite temperatures ? 

It should be noted that a strongly peaked, asymmet- 
ric DOS is a considerably more complex condition for 
ferromagnetism than the Stoner criterion. The latter 
merely asserts that, at T — 0, the critical interaction for 
the instability is determined by the inverse of the DOS 
precisely at the Fermi energy Ep, U c = l/N(Ep). thus 
neglecting antiferromagnetism and the structure of the 
DOS away from Ep. Stoner (i.e., Hartree-FockO) theory 
is a purely static mean-field theory which ignores correla- 
tion effects, e.g., the correlation-induced redistribution of 
momentum states and the dynamic renormalizations of 
the band shape and width. So the question remains: How 
essential are genuine correlation effects for the stability 
of itinerant ferromagnetism at finite temperatures ? 

A third question concerns the suitability of the Hub- 
bard model itself as a model for ferromagnetism. Indeed 
there is no compelling a priori reason why the Hubbard 
model should be a good model for ferromagnetism at 
all. Not only does it neglect band degeneracy, a feature 
observed in all ferromagnetic transition metals (Fe, Co, 
Ni), it also ignores the (weak) direct Heisenberg exchange 
interaction which is equivalent to a ferromagnetic spin- 
spin interaction and, hence favors ferromagnetism in the 
most obvious . aiay.LjEj The proposition by Hirsch and 
coworkersE£t3'tJ that this interaction plays a key role 
in 
al. 



actallic ferromagnetism was disputed by Campbell et 
21 So the controversial question is: 
How important is the direct Heisenberg exchange inter- 
action for the stability of itinerant ferromagnetism in the 
one-band Hubbard model at finite temperatures? 
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In this paper quantitative answers to the three ques- 
tions formulated above are given within the dynamical 
mean field theory (DMFT). The DMFT, a nonperturba- 
tive approach, becppjes exact in the limit of large coordi- 
nation numbers. C3~E3 When applied to d — 3, where the 
coordination number is O(10), the DMFT has proven to 
yield accurate and reliable results_-esEecially in the con- 
text of long-range magnetic order .E3tJ It treats local cor- 
relations exactly while spatial fluctuations are neglected. 
In this situation the momentum integral entering in the 
local propagator will be replaced by an energy integral 
involving only the DOS of the noninteracting electrons. 
The latter may be viewed as an input parameter. In our 
investigation the question concerning the importance of 
the distribution of spectral weight within the band for the 
stability of ferromagnetism will therefore be studied us- 
ing a model DOS of the noninteracting electrons whose 
shape can be changed continuously from symmetric to 
strongly asymmetric by varying an asymmetry parame- 
ter. 

The paper is structured as follows: in Sec. II we present 
the model under investigation, the dynamical mean-field 
equations, and the analytical and numerical steps needed 
to construct magnetic phase diagrams. The model DOS 
is introduced in Sec. III. The results of our investigation 
and quantitative answers to the questions posed above 
are presented in Sec. IV. A discussion where these results 
are put into perspective (Sec. V) closes the presentation. 



II. MODEL AND METHODS OF SOLUTION 

A. Hubbard model with nearest-neighbor exchange 

The minimal model allowing one to treat an asymmet- 
ric DOS, electronic on-site correlations, and the nearest- 
neighbor Heisenberg exchange interaction is given by 

6 = F H ub - 2F Si ■ Sj , where (1) 
-ffHub = - *y (4^„ + h -c-) + U Y ™ijAt (2) 

i,j,a i 

Here S» = | Yl<r<r' ^ia T ^' Cut' w ^ tn tne vec t° r of Pauli 
matrices t. 

We note that there are three other nearest-neighbor 
contributions of the Coulomb interaction which might 
also effect the stability of the ferromagnetic phase (see 
appendix) . 



spin models^ the Hartree-Fock approximation yields the 
exact result for the F-term in high dimensions. 

In the following we investigate the influence of the 
direct exchange term on the properties of the Hub- 
bard model in d — > oo. Since the Hubbard model is 
SU(2) spin-symmetric we can, without loss of general- 
ity, assume a magnetization parallel to the z-axis. The 
HartreeEj decoupling then takes the form Sj • S-,- — + 
(Sf}& + S*(S*) - (Sf){SJ). In terms of the magnetiza- 
tion rh — . rhi /N and its expectation value m = (m) , 

where rhi = 2 S* = hi\ — hii and iV is the number of 
lattice sites, the Hamiltonian (nu can be written asE3 



H = H 



NF* 



NF* 



Hub 



(3) 



Apparently the influence of the exchange term in the 
limit d — > oo is that of a (Weiss) magnetic field which 
vanishes in the paramagnetic phase (m = 0). Therefore, 
in this phase, all one-particle properties of the system are 
those of the pure Hubbard model. However, two-particle 
functions, especial ly th e ferromagnetic susceptibility, are 
modified (see Sec. II C). 

In d = oo, the Hubbard model @ is equivalent to 
an Anderson impurit y . m odel complemented by a self- 
consistency condition£3c3 Written in terms of Matsub- 
ara frequencies oj„, self energy E CT , the DOS of the non- 
interacting electrons N°(e), and thermal average (ipip*)A 
the resulting coupled equations for the Green function in 
the homogeneous phase read: 



G a (iu) n ) = I de 



N°(e) 



(4) 



iuj n + (J, - S ff (icj„) - £ 

— oo 

G a (iuj n ) = -(^ an ^* an ) A . (5) 

The solution of the fc-integrated Dyson equation (^) 
is straightforward and can be performed analytically for 
the DOS used in this paper (see Sec. IV). By contrast 
the solution of Eq. (||) is highly nontrivial (for details 
of the notation see Ref. 0). It is achieved using the 
auxiliary-field quaatum Monte Carlo (QMC) algorithm 
by Hirsch and Fye,E3 where a discretization of imaginary 
time, At = /3/A, is introduced. Here, A denotes the 
number of independent Matsubara frequencies. Physical 
quantities are obtained in the At — > limit. 



C. Calculation of susceptibilities, extrapolation and 
error handling 



B. Dynamical mean- field theory 

Within the DMFT the coupling constants in (|l|) and 
(||) have t o h e scaled with the lattice coordination num- 
ber Z as@I! t = t*/y/Z,F = F*/Z, where we consider 
nearest-neighbor hopping t only. By analogy to classical 



The second order phase transition from a para- to a 
ferromagnetic phase occurs at_±he zero of the inverse 
susceptibility Xf ■> calculatecOo in the paramagnetic 
phase. It is sufficient to perform all simulations for the 
pure Hubbard model since the influence of F* on the 
susceptibility is given by the following random-phase ap- 
proximation (RPA)-like expression: 
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Xf(U,F* 



x f (UA 



(0) 



where \f(U,0, . . .) is the susceptibility of the pure Hub- 
bard model. This type of relation holds for pairs of two- 
particle interactions (here U and F*) in arbitrary dimen- 
sions, whenever one interaction (here F*) is treated in 
HartreCpEbck approximation and the other one (here U) 
exactlyO In general, since the calculation of the suscep- 
tibility involves the derivatrvp-,o£, the self-energy Ej^f* 
with respect to some field this follows from the 

fact that the self-energy of the full Hamiltonian can be 
expressed as 



[Gu,f*] = ^u,o[Gu,f*] + Sq 1 ^. [Gu,F* 



(7) 



Here E[G] refers to the diagrammatic skeleton expan- 
sion of E, where all lines are fully dressed propagators 
G. Since the Hartree-Fock term Eg 1 ]?* only contributes 
in the symmetry-broken phase, all F^-renormalizations 
of Gjj t F* vanish in the symmetric phase. Evaluat- 
ing dSu,F» fdh in the symmetric phase, the first term 
in (Q) leads to the same contributions as without F*- 
interaction, while the second term introduces the RPA- 
like term proportional to F* . For the t-J model in DMFT 
the analogue of @ was derived by Pruschke et al.EJ 

Since F* only enters the calculations via Eq. (0) we are 
left with four physical parameters of the pure Hubbard 
model: Hubbard interaction U, electron density n, tem- 
perature T and an asymmetry parameter a for the kinetic 
energy (see Sec. IV). For each set of these five parame- 
ters Eqs. (|J) and (p|) are iterated with typically 6 x 10 4 
Monte Carlo sweeps until convergence is reached, i.e., the 
difference between two consecutive values of (G a ) _1 = 
(Gcr) -1 — S CT is smaller than 5 x 10~ 4 (measured by the 
norm (2A)' 1 £ CT „ |(^ w ) _1 - (Q%)~ 1 \\ the energy scale 
is defined in Sec. IV). Subsequently eight measurements 
of the susceptibility are performed with a reduced num- 
ber of 2 x 10 4 Monte Carlo sweeps. Thus the result for 
each parameter set consists of an averaged susceptibility 
X/(At) and its statistical error A%/(Ar). We neglect 
the propagation of the error in {Q a )~ x since it is always 
an order of magnitude smaller than A%/(At). The ex- 
trapolation to At = is performed by a quadratic least 
squares fit of x/(Ar), using at least six different values 
of Xf f° r At <E [0.09,0.5]. Further details regarding the 
technical treatment can be found in Refs. |37| , |42] . 

For mean-field theories like the DMFT a linear behav- 
ior of the inverse susceptibility, i.e., a Curie- Weiss law, is 
expected and observed in the vicinity of the transition. 
Thus the Curie temperature Tc can be obtained as the 
zero of a linear fit of X/ 1 ^) drawn from values of Xf 
for four to six different temperatures (see, e.g., Figs. || 
and |, where F* = 2 X$ 1 (§) is plotted). The error of T c 
is obtained from the errors Ax/ (At) by error propaga- 
tion and therefore denotes only statistical, not system- 
atic errors (e.g., due to the extrapolation schemes used). 



However, we checked the accuracy of our results by vary- 
ing the procedure, e.g., extrapolating X/ (At) instead of 
X/(At). 



III. MODEL SPECTRAL FUNCTION 

Due to the vanishing of spatial fluctuations within the 
DMFT the topology of the underlying lattice enters the 
self-consistency equation only via the noninteracting 
DOS, at least for homogeneous phases. The choice of a 
particular model spectral function thus represents a spe- 
cial (not unique) set of hopping elements Uj in the Hamil- 
tonian (0) which characterize the structure of the under- 
lying lattice. Contributions to the kinetic energy by, e.g., 
next-nearest-neighbor hopping can lead to an asymmet- 
rically shaped DOS, which is apparently favorable for 
the stability of ferromagnetism. To investigate this sta- 
bilizing effect quantitatively, we propose a model DOS 
with a shape-controlling parameter. This parameter al- 
lows us to change smoothly from a noninteracting DOS 
with (i) a symmetric shape (mimicking nearest-neighbor 
hopping on a bipartite lattice) to (ii) an asymmetrically 
peaked DOS (similar to a cubic lattice with next-nearest- 
neighbor hopping) to (iii) a DOS with a square-root di- 
vergence at the band edge (e.g., a fee lattice with next- 
nearest- neighbor hopping t' = t/2). The shape of the 
model DOS thus qualitatively captures key features of 
real lattices. 

The spectral function which we use throughout the pa- 
per is given by 



N°(e) 



D 



ae 



(8) 



with c = (1 + VI - a 2 )/(TrD) and bandwidth 2D. The 
well-known semielliptic DOS of the Bethe lattice with 
infinite number of nearest neighbors is recovered for 
a = 0. By increasing the parameter a spectral weight 
is shifted towards the lower band edge (Fig. |l|). For 
a = 1 the DOS diverges at the lower band edge like an 
inverse square- root. The particular choice of the model 
DOS ^ has the advantage that the numerical effort of 
solving the self-consistency equation (Q) is rather small 
since the Hilbert transform can be calculated analyti- 
cally. In the following we set the variance to unity, 
/ de N°(e) e 2 - [J de N°(e) e} 2 = 1, thereby fixing the en- 
ergy scale. This leads to D = 2 for all values of a. For 
a = it is equivalent to choosing t* = 1 on the Bethe 
lattice. 

While for the study of ferromagnetism within the 
DMFT the lattice structure only enters via the DOS it 
is possible to construct (infinitely many) corresponding 
dispersion relations e(fe) or, equivalently, sets of hopping 
elements iy. A realization in d = 1 that is symmetric, 
e(k) = e(—k), and monotonous, de/dk > 0, for e > is 
given by 



3 



a=0 
a=0.5 n 
a=0.9 N°(e) 
a=0.98 
a=1 



N°(e) 



;'\\ 


i i 
a=0.97 




a=0.98 


\\ 


a_n qq 




V — a= 1 








I I 




1.8 

e 



FIG. 1. Model spectral function (|8|) shown for different val- 
ues of the asymmetry parameter a. By increasing a spectral 
weight is shifted towards the lower band edge. The energy 
scale is fixed by setting the variance of the DOS equal to 1. 



/ de'N°(s') =[*L = *® 
J J n it 



(9) 



where an inversion yields s(k). Generalizations to other 
dimensions are possible with, e.g., e(k) — e(\k\)c3 or 

Although in principle one could thus choose a lattice 
corresponding to the DOS (||), this will not be done here, 
since we only study homogeneous phases. Antiferromag- 
netism or incommensurate phases are not expected to be 
important far away from half filling. Only for the case 
a = (Bethe lattice) at n = 1 the stability of an antifer- 
romagnetic phase is investigated. 



IV. RESULTS 

A. The importance of the direct exchange 
interaction 

On a bipartite lattice with perfect nesting, the ground 
state of the pure Hubbard model at half filling is antifer- 
romagnetic for all U > 0, at least in dimensions d > 3. In 
this situation a ferromagnetic state is strongly disfavored 
also in the general model, Eq. ([[]). At large U, how- 
ever, when the model reduces to an effective Heisenberg 
model (which, in high dimensions, is exactly described 
by Weiss mean- field theory), already a small value of 
the direct exchange interaction, F* > 2(i*) 2 /[/,. 
ficient to stabilize a ferromagnetic ground state.l 
Indeed, the Heisenberg model well describes the F*-U 
ground state phase diagram at half filling down to U « 4; 



this is evident from Fig. EL where a comparison with our 
QMC results is shown for a symmetric DOS [Eq. (||) with 
a = 0] . At small U the phase boundary between a param- 
agnetic and a ferromagnetic state is correctly reproduced 
by Hartree-Fock theory (but only for U < 1). This is not 
surprising since in d = oo the F-term is treated exactly 
within this approximation. Also included in Fig. is the 
line below which a fully saturated ferromagnetic state 
becomes unstable, against single spin flips as first com- 
puted by Hirscho for cubic lattices. For the Bethe lat- 
tice this line is given exactly by the Hartree-Fock result 
F* = 4 - U for U < U r = 3 and 



F* — U = — 



(F* +U)(l- y/1 - 16/ (F* + C/) 2 ) 



(10) 



for U > 3. This can be seen from Eqs. 5 and 7 in 
Ref. ^ and the known analytic expression for the Hilbcrt 
transform of the semielliptic DOS. The remarkable agree- 
ment between the QMC results and this curve for U > 3 
suggests that (at zero temperature) the region of par- 
tial polarization is very narrow already at intermediate 
interaction strength U. 

The F*-T phase diagram for a strongly peaked DOS 
(a = 0.98) at filling n = 0.6 is shown in Fig. g. The 
QMC results for the ferromagnetic phase boundary can 
be extrapolated linearly to zero temperature leading to 
a ground state phase diagram. Clearly the values of F* 
necessary to stabilize ferromagnetism are significantly re- 
duced in comparison to the bipartite case. In particular, 
for U — 6 and U = 8 the extrapolation lines cross the 
ordinate at positive temperatures. Thus, an asymmetric 
DOS stabilizes ferromagnetism even in the pure Hubbard 
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FIG. 2. Phase diagram for a symmetric DOS (a — 0): di- 
rect exchange coupling F* vs. Hubbard interaction U at half 
filling (n = 1) extrapolated to T = 0. Open triangles and 
filled circles correspond to the instability of the paramagnetic 
phase (P) against the ferromagnetic (F) and antiferromag- 
netic (AF) order, respectively. Solid line: Hartree-Fock the- 
ory; dashed line: Weiss mean-field theory; dotted line: sinle 
spin flip instability for the saturated ferromagnetic state. 
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model (F* = 0) above a critical interaction strength U c 
with 4 < U c < 6. 

Figs. | and | prove that already a small direct ex- 
change coupling F* can significantly enhance ferromag- 
netic tendencies and thus give the final "kick" towards 
ferromagnetism for systems that are close to an instabil- 
ity. This influence is stronger at larger densities when 
the local magnetic moments are enhanced (Fig. |^). The 
lower critical densities are very small, but larger than 
those predicted by Hartree-Fock theory, since Hartree- 
Fock always overestimates the size of the ferromagnetic 
regime. 



Quite generally, the value of the exchange interac- 
tion F in metals can be expected to be rather small. 
Note, however, that for three-dimensional lattices the 
scaled quantity F* = ZF is an order of magnitude larger 
than the^exchange coupling F itself. Hubbard's crude 
estimated of F/U s» 1/400 would therefore imply that, 
e.g., at U — 6 a (scaled) exchange interaction as large as 
F* = 0.15 is not completely unrealistic. 



B. The importance of the asymmetry of the DOS 
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FIG. 3. F*-T phase diagram for different values of U for a 
strongly peaked DOS (a = 0.98) at a filling n = 0.6. The lin- 
ear extrapolation shows that there exists a critical U c above 
which ferromagnetism is stable even without the direct ex- 
change coupling. 



The dependence of the phase boundary on the asym- 
metry parameter a is systematically studied in Fig. H at 
U = 4 for a relatively small electron density n — 0.3. For 
a symmetric or slightly asymmetric DOS (a = 0, a = 0.5) 
the system only becomes ferromagnetic for F* > 1 even 
at T = 0. For a = 0.9, when the shape of the DOS 
is roughly triangular, the critical F* is considerably re- 
duced. But only when a marked peak develops (i.e., for 
a > 0.95) does the critical F* drop to zero; ferromag- 
netism is then stable even in the pure Hubbard model. 
From now on we restrict our studies to this case (F* = 0) . 

The T vs. n phase diagram is shown in Fig. ^ for 
U = 4 and three different shapes of the DOS ranging 
from strongly asymmetric (a = 0.97) to divergent at the 
lower band edge (a = 1). Evidently the ferromagnetic 
phase is largest at a = 1. We want to stress however, 
that the divergence does not change the physics qualita- 
tively (except at n <C 1). A moderately strong peak near 
the band edge is all that is needed to stabilize ferromag- 
netism. 



0.06 



T 




n 



FIG. 4. T-n phase diagram for different values of the direct 
exchange coupling F* for a strongly peaked DOS (a = 0.98) 
at U = 6. A small direct exchange coupling is seen to en- 
large the stability regime of the ferromagnet, especially for 
densities close to half filling. The lines are guide to the eye 
only. 
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FIG. 5. F*-T phase diagram for different shapes of the 
DOS at U = 4 and n = 0.3. For values of 0.95 < a < 1 
the ferromagnetic phase is stable even without the direct ex- 
change coupling. The lines show a quadratic least-squares fit 
in T. 
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FIG. 6. T-n phase diagram for different shapes of the DOS 
at U = 4. By shifting spectral weight towards the lower band 
edge, the region of stability of the ferromagnetic phase is en- 
larged. The lines are guide to the eye only. 

For symmetric densities of states at U < 12 (Bethe 
lattice and hypercubic lattice) we find suppression of the 
ferromagnetic (as well as the aajtiferromagnetic) suscep- 
tibility away from half filling.E3 This does not exclude 
the possibility for ferromagnetism on bipartite lattices at 
much larger values of U. Indeed, very recently ferromag- 
netism was found within the noncrossing approximation 
for the t—J model on a hypercubic lattice, in the limit 
d — > oo for U > 30 away from half filling.ll3 Apparently, 
at least for moderate U, the bipartite lattice with only 
nearest-neighbor-hopping is not a natural "environment" 
for ferromagnetism - the asymmetry of the noninteract- 
ing DOS is crucial. 



C. The importance of correlations 

In Fig. |^ the T-n phase diagram is shown for different 
values of the on-site interaction U. Evidently the fer- 
romagnetic phase becomes more favorable for increasing 
U: both the maximal Curie temperature and the (upper) 
critical density rise. This effect is seen to be qualitatively 
similar to an increase of the exchange interaction F* or 
the asymmetry of the DOS a (Figs. |^ and^, respectively). 

Our QMC results are compared with Hartree-Fock the- 
ory in Figs. p|-|lO|. We note that, applied to the Hubbard 
model, the DMFT includes Hartree-Fock theory as its 
static limit and is thus superior in any dimension. Fig. || 
shows the vast overestimation of the ferromagnetic phase 
within Hartree-Fock. The maximal Curie temperature 
obtained in this approximation is more than an order of 
magnitude too large. At such high temperatures details 
of the DOS are averaged out and consequently the density 
dependence (e.g., the position of the maximum) is com- 
pletely artificial. In Fig. ^ the U vs. n ground state phase 
diagram is shown and compared to the Stoner criterion.EZI 
At low n the Stoner curve clearly approaches the QMC 
curve. Since the DOS vanishes smoothly at the lower 



band edge for a < 1 both curves diverge for n — > 0. Fig- 
ure [lO] focuses on the limit of large U. The weak coupling 
Hartree-Fock theory fails again: it predicts an unbounded 
linear increase of Tc with U, Tc ~ Un(2 — n)/4, whereas 
QMC shows that Tc has a finite limit for U — > oo. It is 
expected that such a finite limit exists for all densities. 
A saturation is also suggested by the curves in Fig. ^. It 
arises from the suppression of double occupancies by cor- 
relations. In contrast to the Hartree-Fock prediction the 
interaction energy goes to zero for U — > oo, thus only the 
bandwidth remains as an eii&jgy scale. In the special case 
a = 0.98, n = 0.4 one finds^ T C (U = oo) = 0.07 ± 0.02. 




FIG. 7. T~n phase diagram for a strongly peaked DOS 
(a — 0.98) for different values of U. With increasing Hubbard 
interaction the stability regime of the ferromagnetic phase 
becomes larger, especially at higher densities. The lines are 
guide to the eye only. 
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FIG. 8. T-n phase diagram for a strongly peaked DOS 
(a = 0.98) at U = 4. The comparison between Hartree-Fock 
theory (solid line) and DMFT (QMC, circles) reveals the im- 
portance of correlation effects. The dashed line is guide to 
the eye only. 
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FIG. 9. U-n phase diagram at T = for a strongly peaked 
DOS (a = 0.98). The DMFT data (QMC, circles) is ex- 
trapolated from finite temperature calculations. The Stoner 
criterion (solid line) underestimates the critical U c (n) for fer- 
romagnetism, but becomes better at lower densities. The 
dashed line is guide to the eye only. 
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FIG. 10. T-U phase diagram for a strongly peaked DOS 
(a — 0.98) at n = 0.4. The comparison between Hartree-Fock 
theory (solid line) and DMFT (QMC, circles) shows the for- 
mer can describe the Curie temperature Tc(U) neither quan- 
titatively, nor qualitatively. The dashed line is guide to the 
eye only. 

One might argue that comparison of the DMFT re- 
sults should not be made with Hartree-Fock itself but 
with Hartree-Fock plus quantum corrections, since the 
latter are known to reduce many of the deficiencies of 
Hartree-FocL-, Such corrections have been-discussed by 
van DongenE3 and Freericks and JarrelLEil The latter 
authors showed how quantum fluctuations modify the 
Stoner criterion by subtracting the particle-particle sus- 
ceptibility. Evaluating these corrections in the case of 



Fig. |l0| we find that the ferromagnetic phase is com- 
pletely suppressed (as previously observed in Ref. |5l] for 
a symmetric DOS). At a = 0.98 this holds for all densi- 
ties n > 0.3. Thus the "corrections" to Hartree-Fock are 
seen to underestimate the ferromagnetic region by far. 



V. DISCUSSION AND OUTLOOK 

After more than_.three decades of research it has be- 
come clear at lastEjcJ that the Hubbard model can de- 
scribe itinerant ferromagnetism even on regular lattices 
and at moderate [/-values for a wide range of electronic 
densities n. Since ferromagnetism is an intermediate 
to strong coupling problem the question concerning its 
"mechanism" has, in principle, no straightforward an- 
swer. This is in contrast to weak coupling phenomena, 
e.g., conventional superconductivity, which can be ex- 
plained within perturbation theory. Nevertheless a good 
starting point for an understanding of the origin of itin- 
erant ferromagnetism can be obtained in the strong cou- 
pling limit. At U = oo doubly occupied sites are ex- 
cluded and the Hubbard model reduces to a (compli- 
cated) kinetic energy. To avoid doubly occupied sites 
in a paramagnetic phase the DOS is then necessarily 
strongly renormalized compared with the noninteracting 
case, whereas for the saturated ferromagnetic phase the 
interacting DOS is the same as the noninteracting one 
except for a shift between the spin-up and -down bands. 
In this situation details of the structure of the noninter- 
acting DOS become relevant in selecting the state with 
the lowest energy. This physical picture is, in principle, 
similar to that underlying the Nagaoka mechanism; how- 
ever, the latter only addresses the kinetic energy of a 
single hole and it was so far not possible to generalize it 
to thermodynamically relevant densities. Our investiga- 
tions within the DMFT explicitly show that a moderately 
strong peak at the band edge closest to the Fermi energy 
is sufficient to stabilize ferromagnetism. Furthermore a 
strong asymmetry of the DOS implies a nonbipartite lat- 
tice which frustrates the competing antiferromagnetism 
near half filling. 

The mechanism described above is completely differ- 
ent from the mere band shift of the Hartree-Fock the- 
ory. This weak coupling approach does not take into ac- 
count the dynamical renormalization of the DOS in the 
paramagnetic phase and thus predicts ferromagnetism for 
any DOS, even at relatively small values of U and for 
high temperatures. The comparison with DMFT clearly 
shows that Hartree-Fock (i) overestimates transition tem- 
peratures by more than an order of magnitude, (ii) ren- 
ders the dependence of Tc on U qualitatively incorrect, 
and (iii) predicts ferromagnetism for the symmetric DOS, 
where (at least for U < 12) this is not found. These 
shortcomings of Hartree-Fock theory are due to the ne- 
glect of dynamical fluctuations which are at the heart of 
the correlation problem. 



7 



The Heisenberg exchange interaction, not considered 
in the pure Hubbard model, provides another mechanism 
that may order the fluctuating local moments arising by 
the suppression of double occupancies. We found that 
for a symmetric DOS rather large values of F are needed 
to stabilize ferromagnetism. However, for an asymmet- 
ric DOS with a peak near the band edge already small 
values of the exchange interaction may provide the final 
"kick" towards ferromagnetism. In any case it reduces 
the critical on-site interaction and increases the critical 
temperatures of the ferromagnetic phase boundary. 

While the DMFT correctly describes the dynamic fluc- 
tuations of the interacting many-body system, it neglects 
spatial fluctuations and short-range order. Hence one 
should suspect that this approach overestimates the. tran- 
sition temperatures T c . Within DMFT UlmkeEl esti- 
mated Tc for a three-dimensional fee lattice to be of the 
order of 500 to 800 K which is in the range of realistic 
transition temperatures. We may expect spatial fluctua- 
tions to reduce these temperatures. On the other hand, 
band degeneracy, not considered in our model so far, is 
expected to increase Tc- Indeed, band degeneracy and 
Hund's rule couplings which are clearly present in re- 
alistic systems can be rigorously shown to improve the 
stability q f fer romagnetism at least for special parame- 
ter values.ti3ii-3 The incorporation of band degeneracy, for 
which the DMFT also provides a suitable framework, is 
the most important feature that has to he-included in 
future investigations of the Hubbard model.E3 The addi- 
tional nearest-neighbor interactions discussed in the ap- 
pendix may provide yet another mechanism for ferromag- 
netism and will be studied in the future. 
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APPENDIX: NEAREST-NEIGHBOR 
INTERACTIONS 

In Wannier representation the Coulomb interaction 
gives rise to the purely local intera^ipa M as well as 



to four nearest-neighbor interactions ]i!3~E3'L3 Besides the 
Heisenberg exchange interaction these are the density- 
density interaction, the pair- hopphm. term, and the off- 



diagonal "bond-charge-site-charge"lij interaction. The 
latter effectively describes a density-dependent hoppingE2l 
and leads to a narrowing of the band. Hence this term 
is expected to stabilize saturated ferromagnetism. Since 
the quantum dynamics of this term make a systematic 



investigation difficult - even within the DMFT - its de- 
tailed study has to be postponed to the future. The 
pair-hopping.pterm also weakly enhances ferromagnetic 
tendencies. E3E3 

Among all nearest-neighbor interactions the density- 
density term 



H 



NX 



(id) 



(Al) 



is largest and is thus investigated explicitly in the fol- 
lowing. In the case of d-electrons Hubbard roughly es- 
timated this term to V=2-3 eV, an order of magnitude 
smaller than the Hubbard interaction Ua However, since 
there are Z neighbors contributing, the total energy of 
the nearest-neighbor density-density interactions may in 
some materials even surpass that of the Hubbard inter- 
action. This raises the question of the importance of the 
V-term, in particular its influence on the ferromagnetic 
phases investigated in the present paper. ._. 

It was already pointed out by Muller-HartmannEj that 
in the limit d — > oo and with the proper scaling V = 
V*/Z the nearest-neighbor density-density interaction 
reduces to its Hartree contribution which may then be 
viewed as a simple, site-dependent shift of the chemical 
potential. In the absence of broken translational symme- 
try the chemical potential must compensate this shift to 
keep the electron density fixed. Then there is no effect 
at all. 

On bipartite lattices translational symmetry can be 
broken by a charge density wave (CDW) with different 
electron densities on A and B sublattices, i.e., with order 
parameter ticdw — i n A — nB)/2. To study this possible 
ordering we analyze the instability towards a CDW in 
the following. 

Similar to the exchange term F the Hartree contri- 
bution of the interaction V leads, even in the presence 
of other interactions, to an. RPA-like pole in the CDW 
susceptibility [cf. Eq. (|)]£3 



Xcdw(U,V*, . 



Xcdw(U,0, . . .) 
1- V* xcdw(U,0,...) 



Thus a second order phase transition to a CDW occurs 
at V* = 1/xcdw(U,0,...). 

Since next-nearest-neighbor hopping frustrates CDW 
order, the maximal instability towards a CDW is ex- 
pected for the symmetric DOS with a = in Eq. (|§). 
Half filling is optimal in this case. For these parameters 
we determined the phase diagram Fig. 11 employing the 



QMC technique (for details concerning the calculation 
of the CDW susceptibility see Ref. ||). Within DMFT 
a CDW ordering occurs for V* <^ U (at not too high 
temperatures). Compared to the Hartree- Fock approx- 
imation the CDW phase boundary of the full model is 
only slightly moved towards larger values of V* . A sim- 
ilar deviation from Hartree-Fock was-found by means of 
QMC simulations in d — 1 by Hirsch,E3 in d = 2 by Zhang 
and Callaway,E3 and within perturbation theory for both 
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FIG. 11. V*-U phase diagram for the semielliptic DOS 
(a — 0) at T = 0.125. For V* > U a regime with charge den- 
sity wave order (CDW) is established. The antiferromagnetic 
phase vanishes for U > 7.7 ± 0.5 (not shown). The dashed 
line is guide to the eye only. 



weak and st££}M coupling (in arbitrary dimensions) by 



nE'S 



van Dongeni 

All these studies demonstrate that the CDW is stable 
for V* = VZ > U. While this relation may in principle 
hold for some transition metals, ferromagnets apparently 
do not show spatial charge ordering. Therefore the ad- 
equate correlated electron model for a ferromagnet ap- 
pears to have parameters in the range V* < U. Then 
the nearest-neighbor density-density interaction V* has 
no influence on the phase diagram, especially on the bor- 
der of the ferromagnetic phase, at least in d — oo. Even 
in d = 3 the Hartree diagram gives the main contribution 
of the interaction V* since spatial fluctuations, leading to 
genuine correlations, are suppressed as 1/Z . Moreover in 
d = 1 and at half-filling the effect o£J/* on the ferromag- 
netic phase boundary is still small.E3 Therefore, over an 
extended range of parameters the nearest-neighbor term 
V* has almost no influence and thus its importance is 
seen to be much smaller than its value suggests. 
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